Comment on "Infrared freezing of Euclidean QCD observables" 
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Abstract 

Recently, P. M. Brooks and C.J. Maxwell [Phys. Rev. D74 065012 (2006)] claimed that the 
Landau pole of the one-loop coupling at Q 2 = A 2 is absent from the leading one-chain term in a 
skeleton expansion of the Euclidean Adler D function. Moreover, in this approximation one has 
continuity along the Euclidean axis and a smooth infrared freezing, properties known to be satisfied 
by the "true" Adler function. We show that crucial in the derivation of these results is the use of a 
modified Borel summation, which leads simultaneously to the loss of another fundamental property 
of the true Adler function: the analyticity implied by the Kallen-Lehmann representation. 
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I. INTRODUCTION 



In confined gauge theories like QCD, causality and unitarity imply that the Green func- 
tions and the physical amplitudes are analytic functions of the complex energy variables, 
with singularities at the hadronic unitarity thresholds [1]. In particular, the Adler function 
V{Q 2 ) (related to the polarisation amplitude by T>(Q 2 ) = — Q 2 dU/dQ 2 — 1) is a real ana- 
lytic function in the complex Q 2 plane cut along the negative real axis from the threshold 
—AM 2 for hadron production to — oo. This property is implemented by the well-known 
Kallen-Lehmann representation 

AMI 

where TZ(s) is related to the observable cross section cr e + e -^ hajdrons . From ((T]) it follows in 
particular that T>(Q 2 ) is continuous in the Euclidean region Q 2 > and vanishes at Q 2 = 0. 

The renormalization-group improved expansion of the Adler function in massless QCD 
does not satisfy all the properties contained in the above representation. The finite-order 
expansion 

N 

V { S{Q 2 ) = Y. d na n+ \Q 2 ) (2) 

n=0 

is plagued by the unphysical (Landau) pole at Q 2 = A 2 , present in the one loop running 
coupling 

a(Q 2 ) = = - -L— . (3) 
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A modified perturbative QCD series ("analytic perturbation theory"), which imj 
the Kallen-Lehmann representation (pQ) at each finite order, has been proposed in 

Beyond finite orders, the observables can be defined by a summation of the Borel type. 
The Borel transform B(u) of the Adler function has singularities on the real axis of the 
■u-plane Q|: the ultraviolet (UV) renormalons along the range u < — 1, and the infrared 
(IR) renormalons along u > 2 (we adopt the definition of the Borel transform used in Q]). 
While the Borel transform is, for a wide class of functions, uniquely determined once all the 
perturbation expansion coefficients are explicitly given, the determination of the function 
having a given perturbative (asymptotic) expansion is, actually, infinitely ambiguous; not 
only due to the singularities, but because the contour of the Borel-type integral can be also 
varied, without affecting the expansion coefficients of the perturbation series. 
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In Ref. j(J the authors use two different Borel summations of the perturbation series in 
the Euclidean region: for positive coupling, a(Q 2 ) > 0, they choose the integration contour 
along the positive (IR renormalon) axis, 

1 °° 

V PT {Q 2 ) = — [ e -"/W> a « 2 » B{u) du, a(Q 2 ) > 0, (4) 
while for negative coupling the integral is taken instead along the negative (UV renormalon) 



axis: 
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V PT (Q 2 ) = — I e -«/(/3oa(Q 2 )) B{u) dUj a (g2) < Q ^ 



As shown in [6|, the summation based on the above definitions can be expressed as: 

V PT (Q 2 ) = J dru v (r)a(rQ 2 ) , (6) 

5)]. Regulating with the 





in terms of the characteristic function ujd{t) defined by Neubert 
Principal Value the singularity of a(rQ 2 ) at r = A 2 /Q 2 , and taking into account the con- 

n 

tinuity of the characteristic function u>z>(t) at r = 1, the authors of [6( conclude from §6§ 
that the contribution of the leading chain of the skeleton expansion of the Adler function is 
finite and continuous along the whole spacelike axis Q 2 > and approaches a zero limit at 
Q 2 = 0. 

n 

Therefore, in [6| it is shown that by a suitable summation of a class of diagrams in 
perturbative QCD, one recovers a property of the true Adler function, which follows from the 
representation (JTJ). Unfortunately, it turns out that another fundamental property implied 
by same representation (TjQ), namely analyticity in the complex plane, is simultaneously lost. 
In the present Comment, we prove this by calculating the Adler function in the complex 



energy plane with the Borel prescription adopted in 



n 



6|. The calculation uses the technique 



described in [3], based on the inverse Mellin transform of the Borel function. 



II. CHARACTERISTIC FUNCTION AND INVERSE MELLIN TRANSFORM 



As shown in 5j, the function ujx> appearing in ([6]) is the inverse Mellin transform of the 
Borel function B(u): 



Luy 1 

UQ—iOQ 



(7) 



2 



The inverse relation 



B{u) = Jdru v {r) t~ u , 



(8) 



defines the function B(u) in a strip parallel to the imaginary axis with — 1 < Rew < 2, 
where it is assumed to be analytic. 



The function ujt>(t) was calculated in [5| in the large-/?o approximation. The result was 
rederived in [6J. Using (J7]), the calculation is based on residues theorem: for r < 1 the 
integration contour is closed on the right half-u-plane, and the result is the sum over the 
residues of the infrared renormalons; for r > 1 the integration contour is closed on the 
left half-w-plane, and the result contains the residues of the ultraviolet renormalons. The 
reS1 due S of the IR aed UV _ lonS satis* son* Symm e try properties Q, but the, 
contributions are not equal. Therefore uvij) has different analytic expressions, depending 
on whether r is less or greater than 1. Following Ref. jg], we denote the two branches of oox> 
by uj^ and w^, respectively. According to the above discussion, it follows from (JTj) that 



UJ. 



IR, 
V \ 



1 

2vri 



(9) 



du B{u) r u - x - I duB(u) r"- 1 

C+ C- 

where C± are two parallel lines going from to +oo slightly above and below the real positive 
axis, and 
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duB{u)T u ~ l - J duB{u)T u ~ l , (10) 

C' + C'__ 

where C'± are two lines going from to — oo slightly above and below the real negative axis. 
The explicit expressions of and uj^ in the large-/3o approximation are given in Eq. 



(80) of (see also Eq. (2.19) of |7|, where u^ 1 is denoted by and uj^, v by w 
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As shown in [5(, the function ujt>(t) and its first three derivatives are continuous at r = 1. 
Moreover, the explicit expressions given in imply that u^ir) and u>^, v (t) are both 

analytic functions in the r-complex plane cut along the real negative axis r < 0. 



III. ADLER FUNCTION IN THE COMPLEX PLANE 

A closed representation of the Adler function T>pt(Q 2 ) for complex values of Q 2 in terms 
of the characteristic function was derived in [7(. The function T>pt(Q 2 ) was defined for 
large \Q 2 \ by a Borel-Laplace integral along the IR axis, while the expression for low Q 2 
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was obtained by analytical continuation. In the present Comment we use the technique 
presented in 3], adapted for the choice of the Borel-Laplace integral made in 6]. For clarity, 
we shall present the calculation in some detail. 

n 

As in Ref. |6[ we work in the ^-scheme, where all the exponential dependence in the 
Borel-Laplace integrals (jl]) and §5§ is absorbed in the running coupling, and denote by Ay the 
corresponding QCD scale parameter. Let us consider Q 2 complex, first such that \Q 2 \ > Ay. 
Since in this case Rea(Q 2 ) > we use the choice (jl]) of the Borel-Laplace integral with the 
principal value (PV) prescription, taking 

1, 



Vpt(Q 2 ) = ^v (+ \Q 2 ) + v(-\q 2 )], 

where I)W(Q 2 ) are defined as 

p(±)(Q2) = i / e -«/(*«(Q 2 )) B {u) du . 



(11) 



(12) 



Here C± are two parallel lines slightly above and below the real positive axis, introduced 
already in Eq. ([9]). 

Following [7| , we pass from the integrals along the contours C± to integrals along a line 
parallel to the imaginary axis, where the representation (jHJ) is valid. This can be achieved 
by rotating the integration contour from the real to the imaginary axis, provided the con- 
tribution of the circles at infinity is negligible. We consider first a point in the upper half 
of the energy plane, for which Q 2 = \Q 2 \ with a phase < <fi < ir. Taking u = lZe td on a 
large semi-circle of radius TZ, the relevant exponential appearing in the integrals (fl2"|) is 



exp < —TZ 



(\Q 2 \\ 
In , n cos 9 — <b sin 9 
. \ A vJ . 



(13) 



For \Q 2 \ > Ay, the exponential is negligible at large TZ for cos 9 > and sin# < 0, i.e. 
for the fourth quadrant of the complex u-plane. The integration contour defining V^~\Q 2 ) 
can be rotated to the negative imaginary w-axis, where the representation (JSJ) is valid. This 
leads to the double integral 



1 
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du J dr u-d(t) exp 



-u I m — — h it 



A 2 



(14) 
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The order of integrations over r and u can be interchanged, since for positive <fi the integral 



over u is convergent and can be easily performed. Expressed in terms of the complex variable 



Q 2 , the result is 

oo , . oo 

We evaluate now the function V^ + \Q 2 ) given by the integral along the contour C + above 
the real axis. The rotation of the integration contour to the positive imaginary axis is not 
allowed, because along the corresponding quarter of a circle sin# > 0, and the exponent 
( TT3|) does not vanish at infinity for < 0. As explained in [jj, we must perform again a 
rotation to the negative imaginary u axis, for which the contribution of the circle at infinity 
vanishes. But in this rotation the contour crosses the positive real axis, and hence picks 
up the contributions of the IR renormalon singularities located along this line. This can be 
evaluated by comparing the expression (Q of the function uj^-{t) with the definition (j!2p of 
the functions P™: they are connected by the change of variable r = exp[— 1/ (P a(Q 2 ))]. It 
follows that Z>W can be expressed in terms of T>^ > as 

V^=V^ + 2 ^§^(A 2 V /Q 2 ). (16) 

Po Q 2 

The relations (TTT|) . ffTBD and (TIB"]) completely specify the function V PT (Q 2 ) for \Q 2 \ > Ay, 
in the upper half plane Im Q 2 > : 

°° 'A 2 

) PT (Q 2 )= JdrMr)a(rQ 2 ) + j^u;i R (A 2 v /Q 2 ). (17) 
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Using the same method, the function T>pt{Q 2 ) can be calculated in the lower half of the 
energy plane, where Q 2 = \Q 2 \e 1 ^ with — tt < <fi < 0. In this case, the integral along C + 
can be calculated by rotating the contour up to the positive imaginary u axis, while for the 
integration along C_ one must first pass across the real axis and then rotate towards the 
positive imaginary axis. Combining the results, we obtain the following expression for the 
Adler function for complex Q 2 with \Q 2 \ > A 2 ,: 

/ITT A 
dru v (r)a(rQ 2 ) ± - -§ ^ R (A 2 V /Q 2 ) , (18) 
(| Po V 

where the ± signs correspond to ImQ 2 > and ImQ 2 < 0, respectively. We recall that 
the first term in (j!8p is given by the integration with respect to u, while the last term is 
produced by the residues of the infrared renormalons picked up by crossing the positive axis 
of the Borel plane. 



We consider now \Q 2 \ < Ay, when Rea(Q 2 ) < 0. Following j(| we use the definition ([5]) 
of the Borel-Laplace integral along the negative axis. In this case the integral is not defined 
due to the UV renormalons. The Principal Value prescription will be given by fill I) , where 
the are now 

D(±)(Q 2 ) = -L f e -u/(f3 a(Q*)) jB ( M ) dM; ( 19 ) 

Pn J 



C>± 



C± being the two parallel lines above and below the negative w-axis defined in (11 01) . 

We apply then the same techniques as above, by rotating the contours C± towards the 
imaginary axis in the u plane, where the representation (JS]) is valid. If the exponential (TT5j) 
decreases we can make the rotation. If not, we must first cross the real axis and perform 
the rotation. The calculations proceed exactly as before, with the difference that now one 
picks up the contribution of the UV renormalons, according to the relation (fTOj) . This leads 
to the expression of the Adler function for \Q 2 \ < A v 



°" A 2 

V PT (Q 2 ) = ( dru v (r)a(TQ 2 ) ± £ ^^ V (A 2 V /Q 
J p Q 2 



(20) 



where the signs correspond to ImQ 2 > and ImQ 2 < 0, respectively. 

We show now that the limit of the expressions ( fT8l) and (1201) when Q 2 is approaching the 
Euclidean axis coincides with ([6]) . Consider first that Q 2 tends to the real positive axis from 
above, in the region \Q 2 \ > Ay, when T>pt{Q 2 ) has the expression ffl8|) . The integrand has 
a pole at r = Ay/Q 2 . Writing explicitly the real and the imaginary part of the integral we 
obtain, for real Q 2 > A 2 : 



V PT {Q 2 + ie) = Re 



dr ujx>{T)a{TQ 2 



f §[(^(A^ 2 )-^W/Q 2 )]. (21) 



But for Ay/Q 2 < 1, the function uj-d coincides with u^, so the last term in fT2T|) vanishes: 
the imaginary part of the integral in ( fl8l) is exactly compensated by the additional term. 
For Q 2 < Ay, we obtain from (1201) 



VpAQ 2 
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Re 



dr uJx>{r)a{rQ z 



£ % NW) " <W/Q 2 )] , (22) 
Po v 



in the same way. Again the last term in this relation vanishes, since for Ay/Q 2 > 1 the 
function u>d is equal to uj% v . Moreover, one can easily see that the expressions of V PT {Q 2 — 
ie), obtained for Q 2 approaching the Euclidean axis from the lower half plane, differ from 



(I2T|) and (122]) only by the sign in front of the last term, which again vanishes. Thus, for all 
Q 2 > 0, the functions f|T8l) and fl20l) approach the same expression 



V PT (Q 2 ± ie) = Re 



dr uJx>{r)a{rQ' L 



.0 



(23) 



This coincides with the PV regulated integral of the Cauchy type ([6]) which, as shown in Q], 
is finite and satisfies the infrared freezing. Moreover, since ut>(t) is holomorphic (infinitely 
differentiable) for all r > except r = 1, the right-hand side of (l23j) has all derivatives 
defined at Q 2 > 0, except at Q 2 = A 2 , where only the first three derivatives exist This 
means that (1181) and ( 120]) define in fact analytic functions in the regions \Q 2 \ > Ay and 
\Q 2 \ < Ay, respectively. In this way we have obtained, following the approach of Ref. jf|, 
two expressions, (|18p and (120]) . which represent T>pt{Q 2 ) in terms of analytic functions for 
\Q 2 \ > A\r and \Q 2 \ < A\r respectively. 

But the success is illusory, because u^{t) and u^ v (t) are two different analytic functions. 
The expressions (|T8l) and (120]) show that T>pt(Q 2 ) coincides with a certain analytic function 
in the region \Q 2 \ > Ay, but with another analytic function in the region \Q 2 \ < Ay. So, 
the Adler function obtained with the two different Borel representations adopted in 6| is 
not analytic, but only piecewise analytic. This is in evident conflict with the principle of 
analyticity implemented by the Kallen-Lehmann representation ([T|). 

The above discussion refers only to the calculation in perturbation theory. In Ref. 
the authors add a nonperturbative term to the perturbative Adler function. From Eq. (81) 



of 



6[] it follows that the nonperturbative part added to the perturbative function T>pt{Q 2 ) 
given in our relations (fT8l) and (120]) has the form 

V NP {Q 2 )=^uj v {A 2 y/Q 2 ), (24) 

where k is a real constant. Using the fact that cut>(A v /Q 2 ) behaves at small Q 2 
like Q 4 / A v ln(A v / Q 2 ) [5|, one can see from the relations ( 1231) and ( 1211) that the sum 
T>pt(Q 2 ) + T>np(Q 2 ) is finite along the Euclidean axis and vanishes at Q 2 = 0. But it 
fails to be a single analytic function in the complex Q 2 -plane, being only piecewise analytic. 

IV. DISCUSSION 

We have shown by explicit calculation that the Borel prescription adopted in [6] is in 
conflict with analyticity, which is a general property considered fundamental in field theory. 



This result implies that the infrared freezing of the Euclidean observables achieved in [6[] 
has had a price, being possible only at the expense of analyticity. The loss is not only 
of an academic interest: the analytical continuation is the only technique to obtain the 
Minkowskian observables form the Euclidean ones, and all theoretical predictions in field 
theory are based on it. Moreover, the simple model for the complete Adler function proposed 
in 6] cannot represent the physical observable: although it is finite in the Euclidean region 
and exhibits infrared freezing, it is not consistent with the analyticity properties implied by 
the Kallen-Lehmann representation. 

Note also that analyticity is repeatedly invoked by the authors themselves (for instance, 
;he term "analytical continuation" or its verbal analog are mentioned at least eight times in 
6|, in particular in Sections VI and VII, where the Minkowskian ratio 1Z is discussed). Ana- 
lytical continuation is unavoidable even if a smearing procedure is used in the Minkowskian 
region. 

It is worth emphasizing that the result of Ref. |f| is not an intrinsic or natural property 
of the leading one-chain term in the skeleton expansion of QCD, but the consequence of 
a specific, but questionable hypothesis. A step of crucial importance in |6| is the ad-hoc 
redefinition of the Borel integral in the region where the running coupling a(Q 2 ) becomes 
negative. In Ref. 6|, this redefinition originates in a particular utilization of the function 
Ei(z). The authors expressed the Borel integrals, cf. Eqs. (28) and (29) of {(J, in terms 
of Ei(z) depending only on the ratio z = a/z n , where a is the coupling and the z n are the 
positions of renormalons. With the conventional definition of the Principal Value of Ei(z), 
a branch cut is located at a > and z n < 0, or at a < and z n > 0. This implicitly 
selects a specific form of the Borel integral: for a > 0, it is taken along the positive, and 
for a < 0, along the negative real semiaxis, respectively. But this definition is not the only 
possibility. Note that, as pointed out in Ref. p (Section VII), for Q 2 < A v the expression 
(JH]) is not the solution of the renormalization-group equation. We have shown that the use of 
these two different Borel-type integrals defining one single function in two different regions 
is responsible for the loss of analyticity. 

Incidentally, the authors of Ref. |6( admit that the function Ei(z) regulated by the 
Principal Value does not give a reasonable result the for Minkowskian observable 1Z. In 
Section VI they adjust the result by hand, by introducing additional ad-hoc terms (see Eqs. 
(89)-(92) of jd] and 8]). These ambiguous procedures are avoided if analyticity is preserved 



S 



and analytic continuation is performed in a consistent way 
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